MEAN=c(
mean(CD_REAL_1),mean(CD_REAL_2),mean(CD_REAL_3),mean(CD_REAL_4),
mean(CD_REAL_5),mean(CD_REAL_6),mean(CD_REAL_7),mean(CD_REAL_8),
mean(CD_REAL_9),mean(CD_REAL_10),mean(CD_REAL_15),mean(CD_REAL_20),
mean(CD_REAL_25),mean(CD_REAL_30),mean(CD_REAL_35),mean(CD_REAL_40),
mean(CD_REAL_45),mean(CD_REAL_50),mean(CD_REAL_55),mean(CD_REAL_60),
mean(CD_REAL_65),mean(CD_REAL_70),mean(CD_REAL_75),mean(CD_REAL_80),
mean(CD_REAL_85),mean(CD_REAL_90),mean(CD_REAL_95),mean(CD_REAL_100))

UPPER_99=c(
sort(CD_REAL_1)[9950],sort(CD_REAL_2)[9950],sort(CD_REAL_3)[9950],sort(CD_REAL_4)[9950],
sort(CD_REAL_5)[9950],sort(CD_REAL_6)[9950],sort(CD_REAL_7)[9950],sort(CD_REAL_8)[9950],
sort(CD_REAL_9)[9950],sort(CD_REAL_10)[9950],sort(CD_REAL_15)[9950],sort(CD_REAL_20)[9950],
sort(CD_REAL_25)[9950],sort(CD_REAL_30)[9950],sort(CD_REAL_35)[9950],sort(CD_REAL_40)[9950],
sort(CD_REAL_45)[9950],sort(CD_REAL_50)[9950],sort(CD_REAL_55)[9950],sort(CD_REAL_60)[9950],
sort(CD_REAL_65)[9950],sort(CD_REAL_70)[9950],sort(CD_REAL_75)[9950],sort(CD_REAL_80)[9950],
sort(CD_REAL_85)[9950],sort(CD_REAL_90)[9950],sort(CD_REAL_95)[9950],sort(CD_REAL_100)[9950])

LOWER_99=c(
sort(CD_REAL_1)[51],sort(CD_REAL_2)[51],sort(CD_REAL_3)[51],sort(CD_REAL_4)[51],
sort(CD_REAL_5)[51],sort(CD_REAL_6)[51],sort(CD_REAL_7)[51],sort(CD_REAL_8)[51],
sort(CD_REAL_9)[51],sort(CD_REAL_10)[51],sort(CD_REAL_15)[51],sort(CD_REAL_20)[51],
sort(CD_REAL_25)[51],sort(CD_REAL_30)[51],sort(CD_REAL_35)[51],sort(CD_REAL_40)[51],
sort(CD_REAL_45)[51],sort(CD_REAL_50)[51],sort(CD_REAL_55)[51],sort(CD_REAL_60)[51],
sort(CD_REAL_65)[51],sort(CD_REAL_70)[51],sort(CD_REAL_75)[51],sort(CD_REAL_80)[51],
sort(CD_REAL_85)[51],sort(CD_REAL_90)[51],sort(CD_REAL_95)[51],sort(CD_REAL_100)[51])

UPPER_95=c(
sort(CD_REAL_1)[9750],sort(CD_REAL_2)[9750],sort(CD_REAL_3)[9750],sort(CD_REAL_4)[9750],
sort(CD_REAL_5)[9750],sort(CD_REAL_6)[9750],sort(CD_REAL_7)[9750],sort(CD_REAL_8)[9750],
sort(CD_REAL_9)[9750],sort(CD_REAL_10)[9750],sort(CD_REAL_15)[9750],sort(CD_REAL_20)[9750],
sort(CD_REAL_25)[9750],sort(CD_REAL_30)[9750],sort(CD_REAL_35)[9750],sort(CD_REAL_40)[9750],
sort(CD_REAL_45)[9750],sort(CD_REAL_50)[9750],sort(CD_REAL_55)[9750],sort(CD_REAL_60)[9750],
sort(CD_REAL_65)[9750],sort(CD_REAL_70)[9750],sort(CD_REAL_75)[9750],sort(CD_REAL_80)[9750],
sort(CD_REAL_85)[9750],sort(CD_REAL_90)[9750],sort(CD_REAL_95)[9750],sort(CD_REAL_100)[9750])

LOWER_95=c(
sort(CD_REAL_1)[251],sort(CD_REAL_2)[251],sort(CD_REAL_3)[251],sort(CD_REAL_4)[251],
sort(CD_REAL_5)[251],sort(CD_REAL_6)[251],sort(CD_REAL_7)[251],sort(CD_REAL_8)[251],
sort(CD_REAL_9)[251],sort(CD_REAL_10)[251],sort(CD_REAL_15)[251],sort(CD_REAL_20)[251],
sort(CD_REAL_25)[251],sort(CD_REAL_30)[251],sort(CD_REAL_35)[251],sort(CD_REAL_40)[251],
sort(CD_REAL_45)[251],sort(CD_REAL_50)[251],sort(CD_REAL_55)[251],sort(CD_REAL_60)[251],
sort(CD_REAL_65)[251],sort(CD_REAL_70)[251],sort(CD_REAL_75)[251],sort(CD_REAL_80)[251],
sort(CD_REAL_85)[251],sort(CD_REAL_90)[251],sort(CD_REAL_95)[251],sort(CD_REAL_100)[251])

par(mfrow = c(1, 2))

plot(c(1:10, seq(15, 100, 5)), MEAN, type = 'l', lwd = 3, col = 'blue', axes = FALSE, ylim = c(0, .35),
main = '7A. Profile-Based Cliffs Delta for
  1 Through 100 Items per Respondent
  (500 Respondents in Both Groups)',
xlab = 'Number of Items for Both Groups', ylab = 'Profile-Based Cliffs Delta')
axis(side = 1, at = c(1, seq(5, 100, 5))); axis(side = 2, at = seq(0, .35, .05), las=1)
lines(c(1:10, seq(15, 100, 5)), rep(.291068, 28), col='red', lwd=2, lty=2)
legend('bottomright', bty = 'n', c('Reference Cliffs Delta', 'Estimated Cliffs Delta'),
col = c('red', 'blue', 'blue', 'green3'), lty = c(2, 1, 2, 2))

###########################################################################################################

MEAN=c(
mean(PBCD_5),mean(PBCD_10),mean(PBCD_15),mean(PBCD_20),mean(PBCD_25),
mean(PBCD_30),mean(PBCD_35),mean(PBCD_40),mean(PBCD_45),mean(PBCD_50),
mean(PBCD_55),mean(PBCD_60),mean(PBCD_65),mean(PBCD_70),mean(PBCD_75),
mean(PBCD_80),mean(PBCD_85),mean(PBCD_90),mean(PBCD_95),mean(PBCD_100),
mean(PBCD_110),mean(PBCD_120),mean(PBCD_130),mean(PBCD_140),mean(PBCD_150),
mean(PBCD_160),mean(PBCD_170),mean(PBCD_180),mean(PBCD_190),mean(PBCD_200),
mean(PBCD_210),mean(PBCD_220),mean(PBCD_230),mean(PBCD_240),mean(PBCD_250),
mean(PBCD_260),mean(PBCD_270),mean(PBCD_280),mean(PBCD_290),mean(PBCD_300),
mean(PBCD_310),mean(PBCD_320),mean(PBCD_330),mean(PBCD_340),mean(PBCD_350),
mean(PBCD_360),mean(PBCD_370),mean(PBCD_380),mean(PBCD_390),mean(PBCD_400),
mean(PBCD_410),mean(PBCD_420),mean(PBCD_430),mean(PBCD_440),mean(PBCD_450),
mean(PBCD_460),mean(PBCD_470),mean(PBCD_480),mean(PBCD_490),mean(PBCD_500))

UPPER_99=c(
sort(PBCD_5)[9950],sort(PBCD_10)[9950],sort(PBCD_15)[9950],sort(PBCD_20)[9950],sort(PBCD_25)[9950],
sort(PBCD_30)[9950],sort(PBCD_35)[9950],sort(PBCD_40)[9950],sort(PBCD_45)[9950],sort(PBCD_50)[9950],
sort(PBCD_55)[9950],sort(PBCD_60)[9950],sort(PBCD_65)[9950],sort(PBCD_70)[9950],sort(PBCD_75)[9950],
sort(PBCD_80)[9950],sort(PBCD_85)[9950],sort(PBCD_90)[9950],sort(PBCD_95)[9950],sort(PBCD_100)[9950],
sort(PBCD_110)[9950],sort(PBCD_120)[9950],sort(PBCD_130)[9950],sort(PBCD_140)[9950],sort(PBCD_150)[9950],
sort(PBCD_160)[9950],sort(PBCD_170)[9950],sort(PBCD_180)[9950],sort(PBCD_190)[9950],sort(PBCD_200)[9950],
sort(PBCD_210)[9950],sort(PBCD_220)[9950],sort(PBCD_230)[9950],sort(PBCD_240)[9950],sort(PBCD_250)[9950],
sort(PBCD_260)[9950],sort(PBCD_270)[9950],sort(PBCD_280)[9950],sort(PBCD_290)[9950],sort(PBCD_300)[9950],
sort(PBCD_310)[9950],sort(PBCD_320)[9950],sort(PBCD_330)[9950],sort(PBCD_340)[9950],sort(PBCD_350)[9950],
sort(PBCD_360)[9950],sort(PBCD_370)[9950],sort(PBCD_380)[9950],sort(PBCD_390)[9950],sort(PBCD_400)[9950],
sort(PBCD_410)[9950],sort(PBCD_420)[9950],sort(PBCD_430)[9950],sort(PBCD_440)[9950],sort(PBCD_450)[9950],
sort(PBCD_460)[9950],sort(PBCD_470)[9950],sort(PBCD_480)[9950],sort(PBCD_490)[9950],sort(PBCD_500)[9950])

LOWER_99=c(
sort(PBCD_5)[51],sort(PBCD_10)[51],sort(PBCD_15)[51],sort(PBCD_20)[51],sort(PBCD_25)[51],
sort(PBCD_30)[51],sort(PBCD_35)[51],sort(PBCD_40)[51],sort(PBCD_45)[51],sort(PBCD_50)[51],
sort(PBCD_55)[51],sort(PBCD_60)[51],sort(PBCD_65)[51],sort(PBCD_70)[51],sort(PBCD_75)[51],
sort(PBCD_80)[51],sort(PBCD_85)[51],sort(PBCD_90)[51],sort(PBCD_95)[51],sort(PBCD_100)[51],
sort(PBCD_110)[51],sort(PBCD_120)[51],sort(PBCD_130)[51],sort(PBCD_140)[51],sort(PBCD_150)[51],
sort(PBCD_160)[51],sort(PBCD_170)[51],sort(PBCD_180)[51],sort(PBCD_190)[51],sort(PBCD_200)[51],
sort(PBCD_210)[51],sort(PBCD_220)[51],sort(PBCD_230)[51],sort(PBCD_240)[51],sort(PBCD_250)[51],
sort(PBCD_260)[51],sort(PBCD_270)[51],sort(PBCD_280)[51],sort(PBCD_290)[51],sort(PBCD_300)[51],
sort(PBCD_310)[51],sort(PBCD_320)[51],sort(PBCD_330)[51],sort(PBCD_340)[51],sort(PBCD_350)[51],
sort(PBCD_360)[51],sort(PBCD_370)[51],sort(PBCD_380)[51],sort(PBCD_390)[51],sort(PBCD_400)[51],
sort(PBCD_410)[51],sort(PBCD_420)[51],sort(PBCD_430)[51],sort(PBCD_440)[51],sort(PBCD_450)[51],
sort(PBCD_460)[51],sort(PBCD_470)[51],sort(PBCD_480)[51],sort(PBCD_490)[51],sort(PBCD_500)[51])

UPPER_95=c(
sort(PBCD_5)[9750],sort(PBCD_10)[9750],sort(PBCD_15)[9750],sort(PBCD_20)[9750],sort(PBCD_25)[9750],
sort(PBCD_30)[9750],sort(PBCD_35)[9750],sort(PBCD_40)[9750],sort(PBCD_45)[9750],sort(PBCD_50)[9750],
sort(PBCD_55)[9750],sort(PBCD_60)[9750],sort(PBCD_65)[9750],sort(PBCD_70)[9750],sort(PBCD_75)[9750],
sort(PBCD_80)[9750],sort(PBCD_85)[9750],sort(PBCD_90)[9750],sort(PBCD_95)[9750],sort(PBCD_100)[9750],
sort(PBCD_110)[9750],sort(PBCD_120)[9750],sort(PBCD_130)[9750],sort(PBCD_140)[9750],sort(PBCD_150)[9750],
sort(PBCD_160)[9750],sort(PBCD_170)[9750],sort(PBCD_180)[9750],sort(PBCD_190)[9750],sort(PBCD_200)[9750],
sort(PBCD_210)[9750],sort(PBCD_220)[9750],sort(PBCD_230)[9750],sort(PBCD_240)[9750],sort(PBCD_250)[9750],
sort(PBCD_260)[9750],sort(PBCD_270)[9750],sort(PBCD_280)[9750],sort(PBCD_290)[9750],sort(PBCD_300)[9750],
sort(PBCD_310)[9750],sort(PBCD_320)[9750],sort(PBCD_330)[9750],sort(PBCD_340)[9750],sort(PBCD_350)[9750],
sort(PBCD_360)[9750],sort(PBCD_370)[9750],sort(PBCD_380)[9750],sort(PBCD_390)[9750],sort(PBCD_400)[9750],
sort(PBCD_410)[9750],sort(PBCD_420)[9750],sort(PBCD_430)[9750],sort(PBCD_440)[9750],sort(PBCD_450)[9750],
sort(PBCD_460)[9750],sort(PBCD_470)[9750],sort(PBCD_480)[9750],sort(PBCD_490)[9750],sort(PBCD_500)[9750])

LOWER_95=c(
sort(PBCD_5)[251],sort(PBCD_10)[251],sort(PBCD_15)[251],sort(PBCD_20)[251],sort(PBCD_25)[251],
sort(PBCD_30)[251],sort(PBCD_35)[251],sort(PBCD_40)[251],sort(PBCD_45)[251],sort(PBCD_50)[251],
sort(PBCD_55)[251],sort(PBCD_60)[251],sort(PBCD_65)[251],sort(PBCD_70)[251],sort(PBCD_75)[251],
sort(PBCD_80)[251],sort(PBCD_85)[251],sort(PBCD_90)[251],sort(PBCD_95)[251],sort(PBCD_100)[251],
sort(PBCD_110)[251],sort(PBCD_120)[251],sort(PBCD_130)[251],sort(PBCD_140)[251],sort(PBCD_150)[251],
sort(PBCD_160)[251],sort(PBCD_170)[251],sort(PBCD_180)[251],sort(PBCD_190)[251],sort(PBCD_200)[251],
sort(PBCD_210)[251],sort(PBCD_220)[251],sort(PBCD_230)[251],sort(PBCD_240)[251],sort(PBCD_250)[251],
sort(PBCD_260)[251],sort(PBCD_270)[251],sort(PBCD_280)[251],sort(PBCD_290)[251],sort(PBCD_300)[251],
sort(PBCD_310)[251],sort(PBCD_320)[251],sort(PBCD_330)[251],sort(PBCD_340)[251],sort(PBCD_350)[251],
sort(PBCD_360)[251],sort(PBCD_370)[251],sort(PBCD_380)[251],sort(PBCD_390)[251],sort(PBCD_400)[251],
sort(PBCD_410)[251],sort(PBCD_420)[251],sort(PBCD_430)[251],sort(PBCD_440)[251],sort(PBCD_450)[251],
sort(PBCD_460)[251],sort(PBCD_470)[251],sort(PBCD_480)[251],sort(PBCD_490)[251],sort(PBCD_500)[251])

plot(c(seq(5, 100, 5), seq(110, 500, 10)), MEAN, type = 'l', lwd = 3, col = 'blue', axes = FALSE, ylim = c(-1, 1),
main = '7B. Profile-Based Cliffs Delta for
  5 Through 500 Respondents in Both Groups
  (500 Items per Respondent)',
xlab = 'Number of Respondents in Both Groups', ylab = 'Profile-Based Cliffs Delta')
axis(side = 1, at = c(5, seq(50, 500, 50))); axis(side = 2, at = seq(-1, 1, .2), las=1)
lines(c(seq(5, 100, 5), seq(110, 500, 10)), rep(.291068, 60), col='red', lwd=2, lty=2)
lines(c(seq(5, 100, 5), seq(110, 500, 10)), UPPER_99, col = 'blue', lwd=1, lty=2)
lines(c(seq(5, 100, 5), seq(110, 500, 10)), LOWER_99, col = 'blue', lwd=1, lty=2)
lines(c(seq(5, 100, 5), seq(110, 500, 10)), UPPER_95, col = 'green3', lwd=1, lty=2)
lines(c(seq(5, 100, 5), seq(110, 500, 10)), LOWER_95, col = 'green3', lwd=1, lty=2)
lines(c(seq(5, 100, 5), seq(110, 500, 10)), rep(0, 60), col = 'black', lwd=1, lty=2)
legend('bottomright', bty = 'n', c('Reference Cliffs Delta', 'Estimated Cliffs Delta',
'99% Confidence Interval', '95% Confidence Interval'),
col = c('red', 'blue', 'blue', 'green3'), lty = c(2, 1, 2, 2))




